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0^ Abstract 

l/~j In 2+1 dimensions at finite temperature, spontaneous symmetry breaking of global 

^^ symmetries is precluded by large thermal fluctuations of the order parameter. The 

^"^ holographic correspondence implies that analogous effects must also occur in 3+1 di- 

^ ^ mensional theories with gauged symmetries in certain curved spacetimes with horizon. 

p% By performing a one loop computation in the background of a holographic superconduc- 

5—1 
Cd tor, we show that bulk quantum fluctuations wash out the classical order parameter at 

sufficiently large distance scales. The low temperature phase is seen to exhibit algebraic 

long range order. Beyond the specific example we study, holography suggests that IR 

singular quantum fluctuations of the fields and geometry will play an interesting role 

for many 3+1 dimensional asymptotically AdS spacetimes with planar horizon. 



1 Introduction 

The Coleman-Merniin-Wagner-Hohenberg (CMWH) theorem [1] states that continuous 
global symmetries cannot be broken spontaneously in 1+1 dimensions at zero temperature 
or in 2+1 dimensions at finite temperature. The rigorous proofs in these papers implement 
the long appreciated (e.g. [2]) physical fact that divergent infrared thermal or quantum 
fluctuations of the putative order parameter wash out classical expectation values in low 
dimensions. 

The holographic correspondence [3] relates quantum field theories in d+1 dimensions to 
gravitational theories in one higher spatial dimension. The correspondence often admits a 
'large A^' limit in which the dual gravitational theory may be treated semiclassically. In this 
limit the strongly interacting quantum field theory is dually described by fields propagating 
classically on a fixed curved spacetime background in one higher dimension. Each field in 
the 'bulk' corresponds to an operator in the 'boundary' quantum field theory. The bulk 
description of spontaneous symmetry breaking is a field charged under the symmetry which 
has a nonzero purely normalisabla^ behaviour near the boundary [3] . 

The combination of holography and the CMWH theorem implies that there exist space- 
time geometries in one higher dimension than those appearing in the statement of the 
theorem above in which spontaneous symmetry breaking is prevented by strong infrared 
fluctuations. A further difference with the usual CMWH theorem is that the bulk symme- 
try in question is gauged rather than global. The main objective of this paper is to explicitly 
demonstrate this phenomenon via a bulk one loop computation. This will allow us to show 
that the phase of the order parameter is washed out everywhere in the bulk rather than 
only near the asymptotic boundary. It will also enable us to exhibit the expected algebraic 
long range order in the low temperature phase. In general it is difficult to compute loop 
effects in curved spacetimes, where field propagators cannot be found exactly. Our work is 
in the spirit of recent computations |5l El [71 El |9] in 'applied holography' [lOl [11] showing 
how quantum effects due to gapless modes running in loops can be accurately captured. 

A one loop bulk computation is required because in the classical bulk limit the number 
of degrees of freedom 'per site' in the dual field theory becomes very large. In this 'large N^ 



^Near the boundary, the bulk field is characterised by two modes: a 'normalisable' mode that can be 
excited with finite energy and 'non-normahsable' mode that is fixed by the boundary conditions. Given a 
bulk field configuration, the normalisable mode determines the expectation value of the dual operator while 
the presence of a non-normalisable mode indicates that a source for this operator has been added to the 
action of the theory. 



limit, the dangerous fluctuations about the classical 'mean fleld' value are parametrically 
suppressed. This effect was demonstrated explicitly in a field theoretic computation for 
the SU{N) Thirring model in [12j. We will recall more precisely below the sense in which 
the classical limit correctly captures the large N physics, and therefore the sense in which 
previous classical computations using 2+1 dimensional holographic superconductors are 
correct at large but finite A^. While interesting physics resembling the CMWH theorem has 
been reported within a classical gravitational framework [131I14J. the fact that only saddle 
point computations were required in the bulk indicates that those effects were not fluctuation 
driven in the sense that we mean it. Similarly, adding higher derivative '1/A^' terms to the 
bulk action |15] does not move one off a saddle point: local bulk counterterms can capture 
high energy (UV) quantum corrections while the CMWH theorem is a nonanalytic low 
energy (IR) effect. 

Our computation will closely mirror the standard discussions of one loop corrections to 
classical expectation values in flat space quantum field theory. A three point interaction 
in a classical background with broken symmetry will lead to a tadpole diagram describing 
the one loop correction. Sending a hydrodynamic sound mode around the loop produces 
an IR divergent contribution that washes out the phase of the classical order parameter 
and restores the symmetry. This process is illustrated in figure [T] below. Our work is 
complicated by the coupling between several fields and by the need to fix a gauge, but in 
essence is straightforward. 
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Figure 1: The classical background is a planar AdS4 black hole, carrying charge and sur- 
rounded by symmetry breaking scalar 'hair'. We exhibit the divergent tadpole contribution 
of a hydrodynamic sound mode that randomizes the phase of the order parameter. 



We will illustrate the physics of the bulk CMWH theorem by working through a par- 
ticular simple example: holographic superconductors in the probe limit. We note in the 
discussion section, however, that the structure "planar AdS4 horizon + hydrodynamic mode 
= strong IR fluctuations" should be quite general. In particular, there may be cases in which 
the metric itself has strong quantum fluctuations over large spacetime distances at any fi- 
nite value of the AdS radius in Planck units. This would be an interesting and accessible 
quantum gravitational effect that deserves further study. 

2 Classical symmetry breaking solutions 

While much of our discussion will be general, it will be useful to work through a specific 
example. We will consider in detail the case of 'holographic superconductors' with 2-1-1 field 
theory dimensions [161 flTl I18j . A classical treatment of these systems shows spontaneous 
symmetry breaking of a U{1) symmetry at temperatures T < Tc, as we now briefly review. 
We will then go on to show how bulk loop effects will wash out the symmetry breaking 
expectation value over large distances. 

Holographic superconductors are minimally described by a charged scalar field coupled 
to Einstein-Maxwell theory with a negative cosmological constant. For our purposes it 
will be sufficient to work in the probe limit [17J, in which the Abelian-Higgs sector does not 
backreact on the geometry. In this case the background is given by the planar Schwarzschild- 
AdS spacetime 

ds'' = -f{r)dt^ + JM+12 (^^' + ^y') ' (2.1) 



with 



m = i^ (-' - V ) ' (2.2) 



where L is the AdS radius and the horizon radius r+ is related to the temperature by 
T = 3r+/47rL'^. The Abelian-Higgs model then propagates on this background according 
to the action 



S = -^ I d^xy/^ ( -F^ + \d^ - a*|2 + V{\^) ) . (2.3) 



Here e^ is the electromagnetic coupling. The static solutions of interest to these equa- 
tions have a nonzero scalar field ^(r) (taken to be real in the background solution) and 
electrostatic potential <I>(r) = At{r). These satisfy 

r J 

*"+f^ + y')^' = -^^^^ + ^V'{^!). (2.5) 



At this early point for concreteness we will specialise to the specific case 

y = m2|^p. (2.6) 

The form of the potential will not be important in what follows. 



One then looks for solutions to the equations ( 2.4 ) and ( 2.5 ) in which the time component 
of the Maxwell field asymptotes to 

$ = /i + p— + •••. (2.7) 

r 

Here // is the chemical potential which is fixed, while p is the charge density that is deter- 
mined from the equations of motion and regularity at the horizon. It was found in [17] that 
there exists a critical temperature Tc such that for T > Tq the only regular solution has 
vanishing scalar field ^ while for T < Tc the thermodynamically preferred solution has a 
nonzero scalar field with asymptotic behaviour 

^/2 y,3/2+^9/4+L2m2 

This is the 'normalisable' behaviour of a massive scalar field in an asymptotically AdS4 
spacetime and therefore corresponds to a symmetry breaking expectation value (O) of 
the dual charged operator O in the field theory. By numerically solving the background 



equations (2.4) and (2.5) one can obtain e.g. the condensate {0){T) as a function of 
temperature. It was checked that the phase transition at T = Tc is second order with 
mean field exponents. In this paper we will explore the fate of the classical low temperature 
symmetry broken phase under bulk quantum corrections. 

3 Green's functions for multiple fields 

In the remainder of this paper we will proceed to compute a one loop correction to the 
classical bulk expectation value ^(r), whose asymptotic value determines the dual field 



theory expectation value {O) through (2.8). A key ingredient that is needed in any one 
loop computation is of course the propagators of the fiuctuating fields. To this end we now 
obtain the retarded Green's functions for fluctuations about the background. In doing this 
we develop some machinery for computing curved background propagators for flelds that 
mix. 

3.1 Fluctuation equations in Feynman gauge 

There are in general six real flelds to fluctuate: {5At,5Ax^6Ay,5Ar,6'^rT^'^i}^ where the 
last two are the real and imaginary parts of 5^ . The fluctuations are both frequency and 



momentum dependent. Without loss of generality we can take the momentum to be in the 
X direction so that all fields have the form 

<^j{r,x,t) = ^j{r)e-''''+'^\ (3.1) 

Here, as we will often do below, we use <!>/ to denote a generic field fluctuation. It is easily 
seen that the transverse gauge field modes, 6 Ay, decouple from the other modes and so can 
be set to zero for our purposes. 

Our theory has the abelian gauge invariance 

A^A + dX, ^ ^ e'^* . (3.2) 

In order to obtain invertible bulk propagators we need to fix the gauge. It will be convenient 



to use the covariant '^ gauge', which amounts to adding to the action (2.3) the term 

S-gauge = -~ j d'^X./^iVaA'^f . (3.3) 



2e2^ 

With this gauge choice, the ghosts decouple from physical fields and will not play a role 
in our computations. We will immediately specialise to the Feynman gauge ^ = 1 as this 
simplifies the equations of motion. 

From the gauge fixed action we obtain five second order equations of motion for five 
variables (we keep 5 Ay = 0): 



{r-m'.y + (y - ^ - -^ + 'f") m - f^S^r - ^ {r^f^'SAr)' 






a; , , L^k 
-jSAt + — 
J r 

2 r2j,2 



-i^{-6At + ^6A,] =0, 



/ 2 t2i.2 \ 

{r^6A'J - i — ^ - 2^2 J §j^^ + 2* {2^6^r + icoS^i) + iujf'6Ar = , 

// i2 r2/,2 \ / f 



ifSA'J + ( y - ^ - 2^' 1 SA^ + 2ik ( ^5^i + '-6 A, 

^i'^f^Kr + (f - ^ - -H (?)' + /")-•• + »^ (f )' 

-^6At-^^dA, = 0. (3.4) 



These equations were studied in |19j . in the gauge 6Aj. = rather than the Feynman gauge, 
and the hydrodynamic second sound mode was extracted numerically. We shall rederive 
some of their results analytically. While the gauge 6Ar = is often more convenient for 



holographic computations, it also introduces a constraint which complicates loop calcula- 
tions. Mainly for this reason we will stick with the Feynman gauge. 

These equations come together with boundary conditions at the horizon and at the 
boundary. Physical modes should be normalisable at the boundary and ingoing at the 
horizon. It is instructive to characterise the asymptotic behaviours precisely. Near the 
horizon, the general solution takes the form 

SAt -^ at{r - r+Y+\ 5A^ ^ a^^{r - r+y , M^ ~ a^(r - r+)^ , 

6^r^Mr-r+Y , d^ir^ipiir-r+y . (3.5) 

A basis for the five ingoing modes has the following exponents and coefficients 

u = -^: {at,a,„ar,iJr,A} cc {-47rT, 0, 1, 0, 0} , {-21'(r+), 0, 0, 0, 1} , 

{0,1, 0,0,0}, {0,0, 0,1,0}. 
: {at,ax,ar,ipr,ipi} oc {47rr,0, 1, 0,0} . (3.6) 



AttT 

The five outgoing modes at the horizon have the same form but with oj — )■ —uj in the 
above expression for ly and at — )■ —at- We will also use the outgoing modes in constructing 
the Green's function. Recall that 4-kT = ?>rj^/Lp'. Note that the 5A,f = gauge leads to 
behaviours at the horizon that cannot be classified as ingoing or outgoing. This is related 
to the presence of a constraint in this gauge. 

Near the boundary, r — )■ 00, the general solution takes the form 

5 At ~ atr^" , SA^ ~ a^r^ , 5Ar ~ arV^^'^ , 5^^ ~ V'r^'' , 5'^i ~ V"*^^ • (3.7) 

A basis for the five normalisable modes now has 

11 = -I: {at,ax,ar,'ipr,'4^i] oc {eL, 0, 0, 0, 0} , {0, eL, 0, 0, 0} . 
ix = -2: {at,ax,ar,'4)r,'^i] oc {0, 0, eL^,0, 0} . (3.8) 

^ = -3/2 - V9/4 + L'^w? : {at, a^, a,., V'r, V'i} oc {0, 0, 0, eL'^/V2, 0} , {0, 0, 0, 0, eL'^/V2} . 

The orthonormality of these modes with respect to the appropriate inner product will pro- 
vide a simplification shortly, we have chosen the above normalisation for future convenience. 
A basis for the five non-normalisable modes takes the same form, but with the exponents 
changed to 

ft = 0, fi=l, fi = -3/2 + v^9/4 + L'^m? . (3.9) 



The hats are for future use, to distinguish the normahsable from the non-normahsable 
falloffs. For simpHcity we will always work with the quantisation of the scalar field in which 
the normahsable mode is that which falls off faster near the boundary. 

Having identified the physical modes at the horizon and near the boundary, we can 
proceed to write down an expression for the retarded Green's function. 

3.2 The retarded Green's function 

The retarded Green's function satisfies 

KuGfKir,r') = 6iKS{r,r'), (3.10) 

together with ingoing boundary conditions at the horizon and normahsable boundary con- 
ditions at infinityr] In the above expression {/, J, K} run over the five fields and Kjj is the 



second order differential operator of the equations (3.4), written in self-adjoint form, i.e. 



Ku = ^(fu^]+Mjj, (3.11) 



with F,M Hermitian. To put our equations (3.4) in this form we need to redefine one of 
the fields to 

5^i = 6^i- / '^drSAr . (3.12) 



This redefinition will not play a big role in the following, as 6Ar will vanish on the hydro- 
dynamic mode to the order we work. The explicit values of Fjj are given in the appendix. 
In the appendix we derive the following expression for the Green's function of a self- 



adjoint system of the form (3.11) 



, ^f^'^'^r) (W-') ^^j''''^'\r') for r < r' 
G?j(r r') = I ^ \ Jab -^ ^ ^ (3 13) 

<I>j ' [r){W )^fi^j [r ) for r>r 



where the two 'Wronskian' matrices are 

\ dr dr / ' 



The overline denotes complex conjugation. In the above expressions the indices a, b run 
over the ingoing or outgoing solutions at the horizon $('°)'°, $(°"t),a ^^^^ ^.j^g normahsable 



^We will be defining the retarded Green's function by analytic continuation of (the Fourier transform 
of) the Euclidean Green's function Gij{x,y) = —TE{^i{x)^j{y)) from the upper half complex frequency 



plane. In particular this means that there is no i in equation (3.101 



solutions ^(°"y)'°- at the boundary. Note that the complex conjugate of an outgoing solution 
is ingoing (the complex conjugate satisfies the equations when the operator K acts from 
the right). In our particular case, there are five ingoing, five outgoing and five normalisable 



solutions as given in equations (3.6) and (3.8), so the Wronskians are five by five square 



matrices. They are invertible except at special frequencies, as we shall see shortly. The 
Wronskians are constants independent of r, as we show in the appendix. 



The expression (3.13) is a generalisation of the well known form of the Green's function 



for boundary value problems involving a single function. While the result is unsurprising, 
the derivation is not completely trivial and we give details in the appendix. 

Given that the Wronskians are constant, we can evaluate them in the near boundary 
region. To do this we need an expression for the ingoing and outgoing wavefunctions in 
terms of the near boundary wavefunctions. Because the equations are linear, in the near 
boundary region we have 



$ 



(out), a 



ab^{hdy),b 



jjao^ 



-\-V^^ 



ab^{n.n.),b 



$ 



(in), a 



ab^(hdy),b 



jjao^ 



+ v^^ 



ab^{n.n.),b 



+ ■ 



+ ■ 



(3.16) 
(3.17) 



for constant square matrices U, V, U, V. In this expression <I)('^°-).^ are the non-normalisable 

b^Jr^'b and ^^^■''■^'^ ~ ^^r'^o with 
Plugging this expression into 



solutions of (3.9). Recall that near the boundary <I>^ 



(bdy),fe 



coefficients and exponents as described in (3.9) and (3. 



(3.14) gives the leading order behaviour 



W 



ab 



=F^"\Ab 



^6j 



w 



ab 



tV^'"^(A 



^aj 



(3.18) 



The lower sign is for the first mode in (3.8), the rest come with the upper sign. In deriving 



(3.18) we used the fact that we chose our bases such that Fjj(f)'j(f)'^j oc 6"''^, see the explicit 



matrix elements Fjj as given in the appendix, as well as the form of the basis in (3.8). As 
anticipated, the r dependence has cancelled and the Wronskians are constant. 



We can see from (3.17) that whenever there exists a mode that is ingoing at the horizon 



and normalisable at the boundary, then V will be non invertible. It follows from (3.18) that 



the Wronskian is also not invertible. This implies that the Green's function (3.13) will have 



a pole at e.g. the (in general complex) frequency at which the 'on shell' mode appears. 
This is of course the usual statement that Green's functions have poles at the quasinormal 
modes of the black hole, see e.g. [20l[2T] . 

In the following section we proceed to identify the second sound mode and to characterise 
the corresponding pole in the Green's function. 



4 Second sound mode and pole in Green's function 

4.1 Low frequency expansion in the far region and dispersion relation 



In this section we exhibit a gapless solution to the fluctuation equations of motion (3.4). 
That is, we are looking for a solution that exists at small frequencies to and momenta k. 
We expect such a mode to be present as it is the bulk counterpart of the Goldstone mode 
in the boundary field theory. At co = k = we have the exact solution 

S^i = ^, (4.1) 

with all other fields zero. This zero mode is a large gauge transformation of the background. 



The gauge transformation parameter, A in (3.2), is constant and does not vanish at the 



boundary; it therefore induces a global symmetry transformation in the dual field theory. 



The mode (4.1) is however normalisable. It is the normalisable background profile ^ in 
the asymptotically AdS spacetime that allows the gapless sound mode to exist, unlike in 
flat space where the symmetry broken Higgs phase is gapped. At finite but small {co, k} it 
sources other physical modes. We can expand the equations in {u, k} by writing: 

$, = $(0) ^ ^^(l),a; ^ ^^{l),/c ^ ^2^{2W + . . . _ (42) 

This is not a straightforward expansion to perform, because the near horizon (r — t- r+) and 



low frequency {uj — )• 0) limits do not commute. An expansion of the form (4.2) is only valid 
in the 'far' region ^^P 1 ^ ^, away from the black hole horizon. 

(n) 

Within the far region, we now solve for the various $) that are sourced perturbatively 



by the zeroth order solution (4.1), requiring normalisability near the boundary. One finds 



the following correction at linear order in w and k 

SAt = -iu (l - dAl'"""-) , 6A^ = ik(l- 6A]^"'''-) , S^r = icod^r"""- . (4.3) 



The three 'homogeneous' terms are solutions to the co = k = equations. These are chosen 
so that 6 At, 6 Ax and 5^r are normalisable at infinity. The 5 At and 6^r equations are 
coupled: After setting the non-normalisable mode of (^^J?°™' to zero and fixing the overall 
normalisation of 6 At, we are still left with two constants of integration which can be used 
to ensure matching onto the near horizon region and the associated ingoing or outgoing 
boundary conditions. Similarly only the overall scale of 6A]^°'^' must be fixed, leaving a 
constant of integration for matching. We will not need to perform these matchings here. 

(2) 

The next step is to look at the equations for the $]- fields. For our purposes of extracting 
the second sound mode we do not need to solve these equations fully. Rather, we will find 



an obstruction to the existence of the on sheU normahsable mode except on a particular 
dispersion relation (the logic used here is similar to that in |22|). The $]- fields should be 
normahsable at infinity. Using this information together with the form of the first order 



solutions in (4.3), one can expand the equations for the perturbations (3.4) near infinity. 



At leading order one finds that satisfying the 6Ar equation requires 

u^ = elk" , with cl = hm ^"^^^ . (4.4) 

This is the second sound mode found numerically in [19]. That paper also found the 
expected imaginary part of the dispersion relation which is higher order in k. This can be 
extracted in principle from the equations but will not be necessary here. We can also check 



that (4.4) agrees with the (probe limit of the) thermodynamic formula for the second sound 
derived from relativistic hydrodynamics in j23]. Recall that the homogeneous solutions 
are solutions to the lo = k = Q equations whose non-normalisable mode tends to one at 
infinity. The coefficient of the normahsable mode, livur^oo T^^r'^^w™') is then precisly the 
thermodynamic susceptibility for the field theory operator dual to A^u [10]. Thus, in the 
notation of [23J 

Here .^ is a homogeneous background Ax in the dual field theory. The minus sign is due to 
the different sign of the bulk kinetic terms for At and A^ ■ 

In order to find the imaginary part of the dispersion relation we would need to solve 
the equations in the near horizon region, impose ingoing boundary conditions, and match 

in the overlap region ^ <^ ^ 1^1. In our analysis above we noted that there were 

sufficiently many constants to match onto an ingoing solution at the horizon. However, 
we could just as well match onto an outgoing solution. The imaginary part of the sound 
mode pole should have different sign for the retarded and advanced Green's functions (i.e. 
for ingoing and outgoing boundary conditions). The real, non dissipative, part of the low 
frequency dispersion relation should be the same for both Green's functions, consistent with 
what we have found. Via the matching, the ingoing or outoing boundary conditions may 
source other modes in the far region, in addition to those we have considered, which must 
vanish as a; — t- 0. However, because of the linearity of the perturbation equations, these 
extra terms will not interfere with the terms we have found above. 



10 



4.2 Green's function in the vicinity of the pole 

We found in the previous subsection that with either ingoing or outgoing boundary condi- 
tions there is a norniahsable on sheU mode with low frequency dispersion relation uP' = c^k'^. 
Strictly, we showed that this dispersion is necessary to have a normalisable mode at second 
order in frequency and momentum. This is sufficient for the following. We noted that the 
Wronskian matrices will not be invertible at these frequencies. In order to exhibit the pole 
in the inverse Wronskians we need to find solutions away from the zero mode. 



At lowest order in frequency u we have the zero mode (4.1). This is the leading order 



ingoing or outgoing mode and it is also normalisable near the boundary. To obtain the 



Wronskian matrices via (3.18) we need to know the corresponding non-normalisable mode 



generated at higher order in the frequency at the boundary. We noted above that if one 



substitutes the zeroth (4.1) and first (4.3) order solutions into the equations of motion (3.4) 



and assumes that the higher order solutions are normalisable, then the equations can only 
be solved near the boundary when cj^ = c^/c^. If we take w^ 7^ c^k'^, then there is a constant 
term in the 5Ar equation that needs to be cancelled by the addition of a non-normalisable 
falloff for one of the fields at second order. One finds that the non-normalisable mode that 
is sourced is 6^i. Writing as above 5^i = ipir^ we solve for the coefficient to find 

V2d^p 



ipi 



(O) 9/^2 



[oj- 



„2r2 



k^) 



(4.6) 



where we used the expressions in (2.8) and just above (|4.5[) to give the coefficients of 



zero frequency normalisable modes in terms of field theory quantities. The overall scale 



ambiguity of the zero mode in (4.1) will drop out of the Green's function (3.13) 



Now recalling from (3.17) that e.g. $(°^*) ~ U • $('^<^y) + V • ^^f-'^-), we see that the 



matrices V and V, and hence the Wronskians (3.18), have an eigenvalue that is given by 
the right hand side of (4.6). This is because we have found a mode in which both ^j'"^ 



and ^j'^'' in (3.17) are proportional near the boundary (to (5^j),p to leading order in the 



We noted previously that matching to the near horizon region will source additional terms in the far 
region. One might worry that these lead to additional non-normalisable modes. Counting the integration 
constants, the differential equations allow us to set all the non-normalisable terms to zero except for that 
of e.g. S'i'i, and also impose ingoing boundary conditions. At a; = we have an exact solution everywhere 



(4.1 1. At finite co, matching to the near region can only modify the far solution at order cj or higher. In 



particular, <5^i is only modified at order cj or higher. However, we have seen that the order uj" term in S'^i 
in the far region, together with normalisability, completely determined SA^ and SAt at order u, k. These in 
turn determined the non-normalisable term of S'^i at order cu^ ,k'^. Therefore matching to the near region 
cannot change the result ( |4.6[ ). At higher order in frequencies we do expect new non-normalisable terms 
corresponding to the dissipative part of the pole. 
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low frequency limit and close to uj = cj.k . This eigenvalue tends to zero as u — ?• Cgk 



and therefore gives the dominant contribution to the inverse matrix. Plugging into (3.18) 



and then (3.13) and keeping track of the normalisations in (3.6) and (3.8) gives that the 



singular term in the Green's function as w^ — )• c^fc^ is 



Gf....(r,rO = 6^95:^^:1^ + ..., (4.7) 



where the constant 



eHSy'...i^. (« 



a^V 2(2Ao-3) 
Here Aq = 3/2 + y^9/4 + LP"m? is the scaling dimension of the operator O dual to the 



bulk field ^. The expression (4.7) is valid for all r,r' in the far region and does not have 



a discontinuity at r = r' . The numerator in (4.7) will have corrections that go as positive 



powers oi UJ or k. These will not contribute to the long wavelength divergence that we are 
aiming for, as will become clear shortly. Similarly, the other components of the Green's 
function, which we could obtain by turning on different non-normalisable modes, do not 
have a singularity at lJ^ = c^k^ and will also not give IR singular fluctuations. 



The result (4.7) is pretty much the expression one might have guessed from the start. 
Our path has been a little tortuous firstly because we needed to fix the gauge symmetry 
and secondly because we derived an expression for the Green's function for a general set of 



coupled fields. We will outline a shorter, if more formal, route to the answer (4.7) in the 
discussion section. 

5 Washing out the order parameter 



We can now use the Green's function (4.7) to compute the one loop correction to the field 
theory expectation value {O). In fact we can do more than this, we will exhibit a one loop 
IR divergent contribution to the phase of the background field profile ^(r), at all values of 
the radial coordinate r. 

5.1 General formalism for one loop corrections 

Consider the following general Euclidean action in a black hole background with a 3 point 
interaction 

S = J d^xy^) (i$/K/j$j - $/$ja>^A^-^^) . (5.1) 

Here X^'^^ is symmetric in all indices and can depend on the spacetime coordinates. One 
can easily generalise the following to allow for a derivative interaction, but this will not be 
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necessary for us. Expanding the path integral in powers of the interaction, the expectation 
value acquires the one loop tadpole contribution 

{<^Aix)) = ^'}+3J d^yy/g{y)\''^GAi{x, y)GjK{y, y) + ■ ■ ■ ■ (5.2) 

Here Gai{x, y) and so on are the Euclidean Green's function on the black hole background. 
The Green's function will also have a disconnected part, but these do not contribute to the 
IR divergence. We will discuss the regime of validity of the expansion in the interaction 
shortly. Higher point interactions contribute to the tadpole at higher loops and will be 
subleading in the perturbative regime. 

Fourier transform the Green's functions in the field theory directions by writing 

G,j(x, y) = TY,j 0^Gij{r, r', iu^, k)e-'-nXr-r')+iHs-^') (53) 

Here r is imaginary time and the Matsubara thermal frequencies for bosons are 

ujn = 2TTnT , n£Z. (5.4) 



Substituting into (5.2) and looking for a homogeneous expectation value in field theory (i.e. 
with no r or X dependence) we obtain that the one loop correction is 

($A(r)) =tY.I^I dr'y^^)3X''''GAi{r, r' , 0, 0)GjK{r', /, iu^, k) . (5.5) 

The next step is to re-express the Matsubara sum in terms of an integral over real 
frequencies. This is achieved by introducing the contour integral 

T Y, F{iuJn) = ^ [ dz coth ^F{z) . (5.6) 

The contour G has two parts, both running anticlockwise. The first part encircles the poles 
of coth along the upper half imaginary axis (including zero) while the second encircles the 
poles in the lower half plane. The analytic continuation of the Euclidean Green's function 
depends on whether one continues from the upper or lower imaginary axis. The upper 
imaginary axis continues into the retarded Green's function, with singularities in the lower 
half complex frequency plane, while the lower imaginary axis continues into the advanced 
Green's function, with singularities in the upper half plane. This is easy to see from the 
boundary conditions at the black hole horizon [SKUj. Deforming both parts of the contour 
C to the real axis gives 

/oo jQ (' rl'^h r O 

^J ^, J dr'yW)coth^X^'^<GUryAO)ImGUr'y,n,k). 

(5.7) 
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All of the Green's functions appearing here are retarded Green's functions. In deriving 
this expression we used the fact that the symmetric part of the Green's function at real 
frequencies and momenta satisfies 

Gfjj,^ (/, r', n, k) - Gfjj,^ (/, r', 0, A:) = 2ilm Cfj^^ (/, r', Q, k) . (5.8) 

This last equality can be derived from the definition of the retarded Green's function in 



(3.13) using manipulations like those in the appendix. The advanced Green's function G 
is defined via swapping all appearances of <I>(™) and <I)(°"*) in G^. The imaginary part of 
the retarded Green's function is of course (minus) the spectral density, so it is a natural 
quantity to appear here. 

5.2 The IR divergence 



We can now use the final expression (5.7) to evaluate the one loop correction to the classical 
background profile of the scalar field ^ . Using the facts (i) that the singular contribution 
to the Green's function we are interested in appears in the {S'^iS'^i) components as in 



(4.7) and has imaginary part (using the appropriate ie prescription for the retarded Green's 
function) 

Im Gf^^i,^^ (r, r, O, k) = -^e^Q-^{rf sgn{Q)5{n^ - c^k^) + • • • , (5.9) 

(ii) that Sut bj = k = Q the real part of the scalar field 6'^r only mixes with 6 At, as we 



see in the equations of motion (3.4), and (iii) that expanding the action (2.3) leads to the 



following three point coupling involving two 5^jS, 



.5At5^,5^, _ '^^{r) 



^ ^ ^.A..w..w, ^ ^^ ^ (5.10) 



we see that (5.7) implies 



^ 2^y (2^y ^^'V^coth— AG,\,^^(r,r',0,0)ImGf^^,^^(r',r',O,A;) 



27rc 



/d/O?,,,„(r,/,0.0) '-^^ff;^'-' /^..coth^ + 



3e^G f „ , c,k 



■*'^'i^A*~*^^-- '"" 



In the second line we used (5.9) and also introduced a long wavelength cutoff A on the 
momentum integral, which we could think of for instance as a finite volume. As written 
the integral also has short wavelength divergences, but these can be renormalised away and 
are not relevant for the low energy physics we are after. To obtain the third line, we used 

14 



the equation of motion (3.10) for the Green's function Gf^ ^^ , the background equation of 



motion (2.5 ), and integrated by parts twice. A similar argument shows that the electrostatic 
potential does not receive a divergent quantum correction: {5At) = finite. 



The IR singular part of (5.11) behaves as 



{S^r 






(5.12) 



which has a logarithmic divergence as A — )• at all values of the radial coordinate. We 
replaced an arbitrary renormalisation scale in the logarithm with the chemical potential 
fi, which is the natural scale in the system at low temperatures. We can see immediately 
that any terms with extra positive powers of a; or A; appearing in the numerator of the 



Green's function (4.7) will lead to extra powers of k in the integral (5.11) and hence do not 



contribute to the IR singularity. This justifies our neglect of such terms. 



The long wavelength divergence (5.12) is our main result. We should now explain its 



significance and regime of validity. Despite superficial appearances, (5.12) does not imply 



that the correction to the magnitude of the background scalar field diverges, but rather that 
the phase is becomingly completely randomised due to wild quantum fluctuations. If the 
magnitude of the field were diverging, this would be seen in both the real and imaginary 
parts]^ We see above that only the real part of the field has divergent fluctuations. The 
classical expectation value was chosen to be real; therefore as the phase fluctuates around in 
a circle the positive and negative fluctuations of the scalar field in the imaginary direction 
cancel out while there is a net fluctuation in the real part. Alternatively: if the magnitude 
of the field does not diverge, but the real part is becoming ill-defined, it must be because 
the phase is not defined. Thus while the classical solution for the background magnitude of 
the scalar field receives only small quantum corrections, the phase does not have a specific 
value and spontaneous symmetry breaking does not occur. 

For perturbation theory to be applicable, the electromagnetic coupling should be small 



e ^ 1. The one loop correction in (5.12) is down by a factor of e compared to the classical 



value. This remains true so long as the IR momentum cutoff A, which we might think of 
as defining a spatial volume Vol, satisfies 



VolV2 ^ ^-1 ^ ^-^^ ^ ^-l^#/e^ 



(5.13) 



Here # is an order one positive number. Thus, over exponentially large (in e ^ ^ 1) 
volumes of space, in the dual field theory directions {x,y}, perturbation theory is reliable 



*By working with cylindrical coordinates rather than real and imaginary parts, one sees directly that 
there is not a divergent tadpole contribution to the magnitude of the field. 
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and quantum corrections are small. However, on nonperturbatively large volume scales, 
quantum fluctuations eventually dominate and restore the broken symmetry. While the one 
loop correction itself is no longer reliable on these scales, the appearance of the dynamical 
scale Afiuc. should be interpreted analogously to e.g. the perturbative appearance of the 
QCD scale Aqcd^ it is the scale at which a logarithmically running effect becomes order 
one and qualitatively transforms the physics. In our present case it is clear that the physics 
in question is the randomization of the phase. 

In our probe limit, the overall scale of the action has been set by e. However, once 
backreaction onto the metric is incorporated, the magnitude of the on shell action is more 
typically set by some positive power of W, where W measures the degrees of freedom 'per 
site' in the dual field theory. In practice, from a bulk point of view, A^ is some inverse power 
of the AdS radius in Planck units. The limit in which gravitational perturbation theory is 
valid becomes A^ ^ 1, and the lengthscale above which quantum fluctuations dominate the 
geometry becomes 

AfluV = ^-^e#^* • (5.14) 

A potential limitation of the regime of validity of our result comes from the fact that 
we only obtained the Green's function in the far region, ^ <c ^f — 1- Near the horizon, at 



^ < 1 + y, the expression (5.9) cannot be used because the small u) expansion of (4.2) is 



not valid there. However, in the hydrodynamic limit this non-far region gives a vanishingly 



small contribution to the integral over r' in (5.11). This fact combined with the above 



observation that in the far region the integrand is effectively a delta function 6{r — r') 



implies that the answer (5.12) should be true arbitrarily close to the horizon. That is, for 
any given r the dominant values of w and k are sufficiently small that our result holds. 

Before moving on to discuss the presence of algebraic long range order, note a small but 
interesting difference between the computation we have just performed and the more usual 
computation involving a weakly interacting scalar field in 2+1 dimensional flat space with a 



mexican hat type of potential. There, the interaction that leads to a divergence like (5.12) 
comes from the form of the scalar potential. In our case we set the potential to be simply 
a mass term and so there are no such interactions. Instead, we found that the important 
interaction involved the bulk Maxwell field. This Maxwell-scalar interaction is not present 
in the usual field theoretic case, as the symmetry has to be global. This difference highlights 
one sense in which the bulk computation does not fit into the conventional statement of the 
CMWH theorem. 
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5.3 Algebraic long range order 

Our computations so far have shown how, at the level of one point functions, quantum 
fluctuations randomise the phase of the scalar field over large distances and restore the 
classically broken symmetry. Ordered phases are more constructively characterised by the 
existence of 'off-diagonal long range order'. This is the statement that the equal time spatial 
correlator of the order parameter does not vanish at large separations 

Jim {Oix)^O{0)) /O. (5.15) 

\x\—^oo 

In disordered phases, in contrast, the correlator vanishes exponentially at large separations, 
with the correlation length setting the scale beyond which the correlations are small. In 
2+1 dimensional theories to which the CMWH theorem applies, one typically finds that the 
low temperature phase that is classically ordered has instead 'algebraic long range order'. 
We now proceed to show that our holographic system exhibits algebraic long range order, 
both in the bulk and in the dual boundary field theory. Specifically we will find 

lim (O(x)tO(0))~-i^. (5.16) 

Thus in the strict classical limit e — )• 0, genuine long range order is recovered. For any 
finite e the correlation eventually goes to zero at sufficiently large distances. This situation 
is very analogous to that discussed in |12) . At shorter distance scales, which at small e 
can be much larger than /i, the field is locally ordered. In this regime previous results on 
holographic superconductors in 2+1 dimensions will carry over. At larger values of e the 
classical results will be misleading. 



The ingredients going into (5.16) are closely related to those leading to our previous 
result (5.12 ) and so we can be brief. We noted in the previous section that the magnitude of 
the scalar field only receives small quantum corrections, while the phase fluctuates strongly. 
It is therefore useful to write the field in polar form ^ = ^(r)e*^, with ^(r) the classical 
background profile and 9 the phase. The equal time correlator of the order parameter can 
then be written, to leading order in perturbation theory in e, 

(*(r, f)t^(r', 0)) = ^(r)^(r')(e-*(''("'^)-^("''°») = ^(r)^(r')e-5«^("'^^)-'^('-''0))') . (5.17) 

To obtain the second equality one performs a standard manipulation using Wick contrac- 
tions. 



We have already computed the two point correlators of the phase appearing in (5.17). 



For purposes of computing the two point function about a background where the scalar is 

17 



real, fluctuations of the phase are given by fluctuations of the imaginary part of the field 



considered previously: d'^i = ^^(r) 9. Thus from (4.7) the retarded Green's function for the 
phase has the singular behaviour 

GUr,r')= /'® +•••• (5.18) 

Note that this singular term does not depend on r, r' . Using this fact as well as the Fourier 



transformation of (5.3) and the analytic continuation technique of (5.6), the exponent of 



the correlation function (5.17) becomes 

'■oo JO r ^2 



^hj (2^c°th — ImG^,(0,fc)(e*--l) . (5.19) 
In the limit of large separation \x\ — )■ oo, the dominant contribution to the correlation comes 



from the IR singular region of (5.18). To leading order at large separation the integrals above 
can be evaluated to give 

lim mx)em - (0(0)0(0))) = - log^lfl + ... . (5.20) 

We again inserted /i as a characteristic scale. Putting the above together we obtain the 
algebraic long range order 

hm (^(r,f)"^^(r',0)) = ^(r)^(r')|x| ^^ + • • • . (5.21) 

|x|— >oo 

Taking r and r' to the boundary and scaling out the conformal factor, we obtain the 



advertised dual field theory result (5.16). 



6 Discussion 

In this paper we have shown how the quantum fluctuations of a hydrodynamic second sound 
mode in a hairy black hole spacetime wash out the classical symmetry breaking phase. 
This effect was expected as a consequence of the CMWH theorem in the dual field theory. 
Technically our computation depended on identifying the second sound mode, isolating the 
contribution of this mode to the Green's function, and then using standard manipulations 
from finite temperature quantum field theory. These techniques also apply in black hole 
backgrounds due to the boundary conditions at the horizon. 

Before ending with some comments on future directions, we will outline an alternative, 
more direct, route towards obtaining the singular part of the Green's function. This may 
be useful in more general contexts, e.g. with gravitational backreaction, where many modes 
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are coupled. However, the construction is a little formal and should be scrutinised carefully 
before use. The idea, inspired by analogous results for determinants O [2l], is to express 
the Green's function as a sum over quasinormal modes. In terms of the eigenvalues to the 
self-adjoint Euclidean operator 



^ (^Fjj-^^^f^ + MuiiuJn, m'-f^ = Xaiiu^n, ^^f^ , (6.1) 

the Euclidean Green's function may be written as 



Gij{r,r;iu}n,k)=y -^ . (6.2) 

^ Xa{tUJn,k) 

In this expression the eigenfunctions must be taken to be orthonormal. We can now view the 
Green's function as a meromorphic function of complex frequency, analytically continuing 
from u = iujn on the upper imaginary axis to obtain the retarded Green's function. Up to 
'local' UV terms (cf. [24j) the function will be determined by the residues of its poles in the 
complex uj plane. These poles occur at frequencies a;^(A;) at which an eigenvalue vanishes 
X{uJi,{k),k) = 0. These are the quasinormal frequencies of the black hole. At these fre- 
quencies, the eigenfunction ^j {r;iOi,{k),k) becomes the corresponding quasinormal mode 
<&j(r). The retarded Green's function may therefore be written as 

-1 ■^*/'„^^* ('„/^ 



GfAry;u:,k)=J2 "^^ 



MMg.. (6,3) 



If, as in this paper, we are interested in the IR physics due to a hydrodynamic quasinormal 
mode, then we need retain only this term in the above sum over quasinormal modes and 
furthermore in the numerator we can put to = k = 0. In our system, the quasinormal modes 
at w^ = ±Csk together with the eigenvalue 5^A|^^^^ ~ ^k leads to our previously obtained 



form (4.7). Matching the numerical prefactor requires more work, in particular grappling 
with the correct normalisation of the quasinormal modes, and we do not pursue this here. 
An intriguing possibility is that there may be a proliferation of IR singular quantum 
effects in curved spacetime backgrounds in 3-1-1 dimensions with horizons. These need not 
be limited to asymptotically AdS spacetimes or to strictly planar horizons. It may be in- 
teresting to explore cosmological applications |251l26j. The 'membrane paradigm' approach 
to horizons may naturally provide the necessary light hydrodynamic modes. Similarly to 
what we have seen in this paper, sending these modes around quantum loops may lead to 
effects familiar from 2-1-1 dimensional finite temperature physics in fiat space. For instance, 
the shear viscosity diverges in 2-1-1 dimensions [27j . 
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Finally, while we have examined the quantum physics of the low temperature phase, we 
have not looked at the phase transition itself. It seems likely that a bulk renormalisation 
group analysis close to T = Tq will show a Berezinskii-Kosterlitz-Thouless (BKT) transition 
|28j rather than the classical second order phase transition. It would be nice to exhibit 
various BKT scalings explicitly. 
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A Green's function for one dimensional boundary value prob- 
lem with multiple fields 



In this appendix we derive the formula (3.13) for the Green's function as well as the state- 



ment that the 'Wronskians' (3.14) are constant. That is, we will obtain the Green's function 
for a one dimensional boundary value problem involving multiple fields. 
The Green's function satisfies 



Fij^\ + Mu\ G'iKir, r') = 5iK6{r, r') , (A.l) 



d f ^ d 
dr 



with F, M Hermitian matrices, together with boundary conditions. To tie in with the main 
text, we will call the boundary conditions at r = r+, 'ingoing/outgoing', and the boundary 
condition at r — )• c«, 'normalisable'. Also, we will assume that if there are A^ fields, i.e. 
/ = 1 . . . A^, then there are A^ ingoing boundary modes and A^ normalisable modes. We will 
also need the A^ outgoing modes whose complex conjugates are ingoing. We can therefore 
define the N x N square matrices 

<&Wout = ^If °"*\ ^oo = <l>ff^\ G^^Gfj, P^Pa,, Q^Qat. (A.2) 

Here the matrices ^^^j"^ and ^^^j'^ combine the A^ field components of the A^ ingoing, 
outgoing and normalisable solutions, each labelled by a = 1 . . . N. We then make the 
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following ansatz for the Green's function 



G^{r,r') 



^Ini^)Pi'''')^oo{r') for r < r' 
^'Lir)Q{r')^ont{r') for r > r' 



(A.3) 



We now proceed to solve for P, Q and to show that they are in fact constant and independent 
of r' . As in the text, over line denotes complex conjugation. 

There are two matching conditions to impose at r = r' . The fields must be continuous 



and the derivative of the fields must jump according to (A.l). Continuity is immediately 
seen to imply 

P = (<l>oo'J>rn')^Q^out^;o' • (A.4) 

The discontinuity of the derivative at r = r' then implies, using the previous continuity 
equation to eliminate P, 






in^-lT^T 



\ dr 
Let us introduce the two matrices 



dr 



^^^''^L]Q = u. 



(A.5) 

(A.6) 
(A.7) 

(A.8) 
The advantage of this expression is that we will now show that W is constant and that S 



S 
W 



dr 
dr 



dr "^' 
dr °°' 



and thereby rewrite (A.5) as 



-5$r^i^C + ^)Q = Id. 



vanishes. Given that W is the first of the Wronskians as defined in (3.14) in the main text, 



this gives the result we are seeking, Q = W ^. A few steps of algebra, using the vanishing 
of 



s = ^^f"'- ""- 



dr 



dr ~' 



then show that P = W ^. Here 



— d^c>o T — d^T„ 



dr 



dr 



(A.9) 



(A.IO) 



completing the derivation of (3.13). 



Constancy of W, W and S, S is immediate from acting with ^ and using the Hermiticity 
of F and M. Given that S and S are constant, to show that they vanish it is sufficient 
to show that they vanish at any point. The natural points to take are the horizon for 
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S and infinity for S. For S, vanishing follows immediately from the orthogonality of the 



normalisable modes in (3.8). For S we must work a little harder. Firstly, taking into account 



the shift (3.12), we have 



e'Frr = ^ + ?4^ (/; ^dr)' , e^F 5^ = ?^ /^ ^dr . (A.ll) 



Using the ingoing modes of (3.6) we can compute S on each pair of modes. It is found 
that all components of S vanish at r = r+. This vanishing depended on the combination 
of <1>?^ and $out appearing in our ansatz (A.3). The combination can be deduced from first 



principles by working in Euclidean signature, where there is a Dirichlet boundary condition 
at the horizon, and then analytically continuing to obtain the retarded Green's function. 

D 
The vanishing of S and S is presumably general for self-adjoint boundary value problems, 
as it is necessary for the Green's function to have the required symmetry properties under 
r -^ r' . For Dirichlet boundary conditions, for instance, supposing the Wronskians are 
constant and finite, then the S and S matrices must be constant and zero because they 
have a faster falloff near the relevant boundary. 
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